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Models with higher order derivative terms in the kinetic energy appear not only as effective 
theories, they can be considered as elementary, renormalizable models in their own right. The 
extension of Higgs mechanism is discussed for Yang-Mills-Higgs models where the kinetic energy of 
the scalar field contains higher order derivatives. The connection term of the covariant derivatives 
generates a local potential for the gauge field in the phase with spontaneously broken global gauge 
symmetry which may break Lorentz symmetry in the vacuum. Relativistic, massless gauge bosons 
spectrum is found for small fluctuations where the only remnant of the breakdown of the Lorentz 
symmetry is a shift of the energy induced by the expectation value of the temporal component of 
the gauge field, reminiscent of a chemical potential. 
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There are two different ways in which higher order derivative terms may arise in the Lagrangian of field theoretical 
models. The usual way is to start with an elementary, renormalizable theory and look for the effective dynamics, below 
some heavy mass particle production threshold. The effective Lagrangian usually contains space-time derivatives of 
Q-i arbitrary orders, reflecting the nonlocality of the dynamics at the characteristic scale of the heavy particle. Another 
possibility is to include some higher order derivatives in the Lagrangian of an elementary, renormalizable model as a 
part of the dynamics or simply as a regulator [lJ-Q ■ The development of both kinds of models is seriously hindered 
by notorious problems with unitarity [8[ . Though this is a trivial problem in the case of the effective theory, arising 
from the incompatibility of a smooth UV cutoff in energy and the truncation of the gradient expansion, there is no 
satisfactory solution in sig ht @. 

The imposition of the most natural condition to ensure unitarity, namely to keep the elementary excitation spectrum 
real, leads to an unexpected difficulty: theories with higher order derivatives handle states with positive, negative and 
vanishing norm. States with negative norm defy physical interpretation, the weight of states with vanishing norm 
remains unbounded even for unitary time dependence and may generate instabilities. But the unitarity of the time 
evolution can be confirmed in PT-invariant theories because the signature of the norm is given by the PT parity 

mm- 

Once the troubling aspects of the higher order derivative terms can be controlled it is natural to inquire if such 
theories can display new type of phenomenology. There are two features of theories with higher order derivatives 
which make them unique. One of them is a very economical unification. It is well known in solid state physics 
that higher order time derivatives in the kinetic energy lead to several bands. In the present context of relativistic 
models the higher order space-time derivatives in the kinetic term generate several relativistic particles, preons. Two 
nontrivial manifestations of preons in the presence of spontaneously broken continuous symmetry are considered in 
this work. It is shown below that the spontaneous breakdown of a global, continuous symmetry leads to massless 
Goldstone preons as expected, leaving the possibility for the remaining preons of the Goldstone mode to be massive. 
In case of a local symmetry the Higgs mechanism is modified. The gauge field remains massless and the usual balance 
between the disappearing Goldstone modes and the appearing longitudinal modes of the gauge field, observed in the 
usual Higgs mechanism does not hold anymore. In particular, the number of Goldstone preons of the scalar field 
which are removed from the physical, gauge invariant sector of the theory is not necessarily related to the number of 
preons appearing in the massless gauge field. 

Another unusual possibility, offered by the higher order covariant derivative terms in gauge theories is the sponta- 
neous breakdown of Lorentz symmetry. The status of Lorentz symmetry is remarkably nontrivial in particle physics. 
Einstein's special relativity relies on macroscopic objects like meter rods and watches and is lacking its solid observa- 
tional basis in the microscopic, quantum regime. The rules of formal implementation of Lorentz symmetry in QED 
are confirmed by the high precision measurements but their relation to Einstein's special relativity is less clear in 
the absence of measuring devices of velocity of propagation within an atom or a nucleon. The retarded propaga- 
tors, describing the linear response to an external source, are supposed to satisfy causality and are made to vanish for 
space-like separation. But the Feynman propagator is nonvanishing for space-like separation allowing acausal coupling 
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among quantum fluctuations, a possible source of acausal higher order, nonlinear response. 

This issue can be raised in a more phenomenological manner, as well, by recalling an important lesson of the 
renormalization group method: the measured values depend on the scale of observation. In other words there are 
no constants in physics. Any experimental constant is actually a relatively large plateau with slowly varying scale 
dependence only. Is the velocity of light an exception to this rule? Lorcntz symmetry is obviously lost at long 
distances due to the expansion of the Universe. Since the experimental verification of special relativity is difficult 
at extremely short distances one is tempted to look for possible quantum effects, generating different speed of light 
at microscopic scales without unsettling Einstein's special relativity stated and experimentally verified beyond any 
doubt at intermediate scales. 

Actually, one stumbles on a problem with Lorentz symmetry at short distances independently of the suspicion raised 
by renormalization group. It is the impossibility of constructing UV regulator with relativistic symmetries. The reason 
is the impossibility of defining a region in the energy- momentum space in terms of invariant length due to the divergent 
volume of the Lorentz group. The regulators, based on Wick-rotation always contain an O (1/A) Lorentz violating part, 
A being the energy scale of the UV cutoff. This is not a problem for renormalizable non-gravitational theories where 
the global Lorentz symmetry is recovered when the cutoff is removed [l2l ] but poses a serious problem for quantum 
gravity where the gauge symmetry, involving Lorentz transformations must be implemented in a cutoff-independent 
manner. The way out might be either the use of stringy description as a regulator or giving up altogether the Lorentz 



structure at the cutoff [13L_ 14L The recent observations about the possibility of generalizing special relativity to 



include the Planck length |15l Il6l | represent a promising step in constructing an interpolating theory. 

Therefore one turns to the question of the breakdown and the restoration of Lorentz symmetry in the function 
of the scale of observation in quantum field theory. The simplest way to break a symmetry is to introduce explicit, 
non-symmetrical terms in the action. But it may happen that a symmetry, present formally on the level of the bare 
action is broken dynamically at some finite scale. If the symmetry breakin g is driven by asymptotically long distance 
modes then it is called spontaneous. It was found a long time ago fl7l Il8|| that the explicit breaking of Lorentz 
invariance, introduced by the specific choice of metric tensor or the vierbein at the cutoff scale, is suppressed at finite, 
physical scales as the cutoff is removed. If the Lorentz symmetry is broken explicitly in a Yukawa model at the cutoff 
scale by non-derivative couplings then it is found to be restored at finite energies after the cutoff is removed [19|. 
A particular symmetry breaking, characterized by a four-vector in a Chern-Simons action has been proposed in a 
CPT-even manner [20| . Another symmetry breaking pattern, achieved by dimension-5 operators has been identified in 
Ref. pH . The issue of generating Lorentz breaking terms by radiative corrections from simpler theory with CPT-odd 
and non-Lorentz- invariant terms in its action has been considered as well [22j - [25| . 

The astrophysical long distances violation of Lorentz symmetry hints at a spontaneous s ymm etry breaking, a 
gravitational Higgs mechanis m [2 (1 and relations between gravitons and Goldstone modes (27l. |28|. a scenario based 
on an earlier, Abelian idea [29h31| |. One can imagine non-gravitational roots of the spontaneous breakdown of Lorentz 
symmetry. In a generic relativistic model the nonvanishing vacuum expectation value of a vector or tensor field signals 
such a symmetry breaking (2(| |32[ . This order parameter is not necessarily an elementary field, it can be generated 
dynamically (33j , as well. It should be noted that the spontaneous breakdown of Lorentz symmetry might by hidden 
from observation, a point of view which leads to the systematic establishment of gauge theories [34j . 

A simple dynamical breakdown of Lorentz symmetry has already been demonstrated in an effective scalar QED fTll j . 
When higher order covariant derivatives occur in the kinetic energy of a charged scalar particle then the connection 
term of the covariant derivative induces a local potential for the gauge field and this may lead to nonvanishing 
vacuum expectation value for the gauge field. The excitation spectrum is not necessarily relativistically invariant in 
this phase therefore one has to distinguish carefully the relativistic, massless preon spectrum from the relativistically 
noninvariant, gapless spectrum. In case of scalar QED the gauge field is the Goldstone mode, arising from the 
spontaneous breakdown of the symmetry with respect to Lorentz boosts and remains gapless. The usual, relativistic 
Green functions are recovered for the gauge field, pushing the signature of the breakdown of the Lorentz symmetry 
beyond the tree-level dynamics of the photon field. We consider a non- Abelian Yang-Mills-Higgs model in this work 
where the naive argument would allow the mass generation for the gauge field in the generalized Higgs mechanism. 
Instead, the tree-level dynamics of the gauge field remains gapless and the only breakdown of Lorentz invariance is 
the emergence of a shift of the energy of the charged gauge field components, reminiscent of a chemical potential, 
provided by the expectation value of the temporal component of the gauge field. 

The local potential, generated by the connection term of the covariant derivatives in the spontaneously broken 



phase of the global gauge symmetry is reminiscent of the bumblebee models 32[. The Lorentz symmetry is then 
recovered on the tree [35| and one- loop level [3(1 H3] and with non- Abelian symmetry (38| . This was argued to be a 
mechanism to render the breakdown of Lorentz symmetry invisible {39] ■ Our model differs from those used in these 
considerations because our local potential for the gauge field leaves the invariance under local gauge transformations 
intact. 

The preon content of a scalar field with a global internal symmetry is discussed in Section [TTl starting with the well 
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known difficulties, caused by the higher order derivatives in the kinetic energy in Section III Al The preon fields and 
condensate are introduced in Section III Bl and III CI respectively. The extension of Goldstone theorem for preons is 
presented in the simple context of an effective potential in Section lllDI The extension to local symmetry is addressed in 
Section ITO1 where first the tree-level vacuum is identified in Section [ill Al followed by the calculation of the quadratic 
action for small fluctuations in Section llll Bl The gapless nature of the scalar sector is verified in Section llll CI Section 
IIVI contains a more detailed discussion of a simpler case where the gauge group is chosen to be U(2) = SU(2) x U(l). 
The gauge field condensate can be either purely non-Abelian or completely Abelian in this case, these possibilities are 
considered in Sections IIV Al and IIV Bl respectively. Section IVl contains our conclusions. The details of the calculation 
of the quadratic action are presented in the Appendix. 



II. PARTICLES IN THEORIES WITH HIGHER ORDER DERIVATIVES 



The higher order derivative terms in the Lagrangian of a quantum field theory describe less local coupling at the 
level of the UV cutoff. This can be seen easily in lattice regularization where a finite difference of order n in the action 
represents coupling of field variables, placed n times the lattice spacing apart and renders the dynamics nonlocal up 
to these distances. Such a short distance modification of the dynamics used to arise in two different manners. One 
way is to implement a regulator which is kept as a physical law rather than being removed in the renormalization 
process. In fact, higher order derivatives in the kinetic energy have been used as possible regulators for a long time 
[lHa] . Another circumstance is when one is interested in the physics at energies well below the energy of the opening 
of particle creation channels and seeks an effective theory by eliminating the particle in question. The perturbative 
construction of such effective theories yields effective bare action with non-polynomial dependence in the space-time 
derivatives, a class of theories beyond our analytic capabilities. The truncation of the dependence on the space-time 
derivatives of the action into a polynomial of finite order renders the theory soluble within the scheme of perturbation 
expansion. 

The problem of theories with polynomial higher order derivatives, addressed in this Section, is a grave conflict with 
the probability interpretation of quantum mechanics. 



A. Difficulties 



The source of the problems in theories with higher order derivatives is the unusually large set of frequencies which 
solve the linearized equation of motion. Let us consider this problem in the framework of a complex scalar field 4>j(x), 
j = 1, . . . , n whose dynamics is defined by the Lagrangian 

L = ^L(-n)<p-v(ct>U), (i) 

characterized by two real functions, the local potential V(<ft<j>) an d L(p 2 ). The interaction, V"((ft4>), where the 
prime denotes the derivative with respect to the argument, (ft '</>, is supposed to be sufficiently weak to justify the 
discussion below and the kinetic energy, L(p 2 ), is assumed to be a polynomial of order rid- The mass term is placed 
into the potential by requiring L(0) = 0. Let us assume a homogeneous external source, coupled linearly to the field 
which sets 4>(x) = 4>. The dispersion relation of infinitesimal fluctuations around this state is given by the equation 
det[L(p 2 ) -M 2 } = with 

M 2 = tv'($4>) + 4> ® ipv"(fi4). (2) 

Owing to the higher order derivative terms in the kinetic energy we have 2nd different values of the frequency, 
corresponding to a fixed three-momentum p. The generalization of usual construction of the Hamiltonian in classical 
mechanics [40( leads to nondefinite kinetic energy in the Hamiltonian. 

The boundedness of the Hamiltonian from below can be saved in quantum mechanics by assuming negative norm 
in the subspace where the negative contributions of the kinetic energy are acting [8]. This can be seen easily in the 
case of the Lagrangian ([1]), whose propagator, 



D ~L(p 2 )-M 2 (3) 



is given by the partial fraction decomposition as 



L(p 2 ) -M 2 ^ p 2 



^ " IT, (4) 
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in terms of the spectrum {m^}, the set of the roots of det[L(p 2 ) — M 2 }, and z n = [dL(p 2 )/dp 2 ]~ 1 \ p 2 =m ^. It was 
assumed in deriving Eq.(j4|) that each eigenvalue A a (p 2 ) of L(p 2 ) — M 2 has single roots only. If the root p 2 = m 2 of an 
eigenvalue is of £-th order then the right hand side may contain terms Zj^/(p 2 — m 2 ) k with 1 < k < £. It is attractive 
to interpret this form of the propagator by assuming that the field <fr handles the propagation of several particles. 
Though the constants z n are real, they can be negative. In fact, any continuously derivable, real function, e.g. 
x — > A a (x), displays alternating sign in its derivative as it crosses zero at odd order roots. If the negative contribution 
of "particles" with z n < to the propagation is interpreted as the result of the negative norm of their state then their 
contribution to the kinetic energy becomes positive. Another problem, related to the unusual frequency spectrum is 
that poles may be complex and the unitarity of the time evolution is lost in this case. 

One can argue that the negative norm states and the nonunitary amplitudes could be suppressed by a sufficiently 
low UV energy cutoff [7|. But once negative norm states appear in the theory there are also states of vanishing 
norm and their virtual contribution to the transition amplitude diverges exponentially in time, free of any observable 
constraint. These modes can be eliminated by imposing appropriate boundary conditions at the final time Q but 
they in turn generate acausality and nonunitary time evolution, closing the vicious circle. The perturbative attempt 
to restore unitarity by modifying Feynman rules @ runs into difficult ies, as well. 

Negative norm states have already been used in gauge theories 0, H2| , but these states can be excluded from the 
final state, as long as the system starts in the physically interpretable positive norm subspace, by means of gauge 
invariance. A similar symmetry argument was found useful to reach similar conclusion covering matter fields as 
well. One can show that the signature of norm of an excitation, created by the elementary physical fields is identical 
with product of the time and space reversal parities of the state [HI HH ■ Hence the unitarity of the time evolution 
follows within the positive norm subspace in time reversal invariant dynamics as long as the values rr? n are real. It is 
interesting to note that this argument is valid in quantum mechanics only. In classical physics a discrete symmetry 
leads to conservation law for harmonic systems only where the equation of motion is linear. 

Higher order derivatives always occur in the effective bare action when the dynamics is sought below a particle 
creation threshold and their leading problems, mentioned above arise as soon as the effective action is truncated to 
polynomials of finite order in the gradient expansion. This is clearly a problem caused by the approximation only, 
the imposition of a UV cutoff below the heavy particle threshold solves all problems when the full effective action is 
used. But the systematic treatment of a sharp energy cutoff is rather difficult and one is usually satisfied by treating 
the higher order derivative terms as perturbation within the gradient expansion. 



B. Preons 



The identification of a solution of a classical equation of motion, a differential equation of the order 2nd in the 
time derivatives needs 2nd parameters. This suggests that the field <f> of the theory (fT]) describes particles, called 
pre-particles or preons below. The preon level description of this theory is defined by the partition function 

"* n 

where it is the sum of the preon fields which enters into the vertices and which is coupled to the external source. This 
sum, denoted by <f>, can be introduced by a constraint, represented by the insertion of 1 into the path integral, 

e'^^l = f D^Mn^^WnW 1 ) _ ^^ n ( a .)) e */«faE»to l *n-vwV)+i^+iV] < ( 6 ) 

n x n 

The Fourier integral representation of the Dirac-delta yields the form 
jwvj] = f D [ a ] D [ a 1] D [tf] D [flY[D[<li]D[^ (7) 



e 



where the preon field can easily be integrated over, 

e'^i'l = J J}[ a ]fl[ a t]J)[^j)[^'/ i [-« t E,»-.«^W t *)+(« t +] t )W t («+ ) )] ] (8 ) 

and one finds 

e iW[j,i + ] = f D ^} D ^ e i f dxi^D-^-V^UHjU+^j} ^ (9) 
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with 

D = J2D n , (io) 

n 

after eliminating the auxiliary Fourier variables. 

The poles of the propagator in energy at vanishing momentum define the preon mass spectrum and we have rid 
preons if the function L(p 2 ) has single roots only. The inequalities m 2 > and z n > are needed to guarantee unitary 
time evolution and positive definite norm for preons, respectively. In case of higher order roots the corresponding 
preon gives no contribution to scattering amplitude, it is "confined" . 

The preons are not elementary constituents of the excitation of the field 0, the latter can rather be imagined as 
appearing in the form of one of the preons, the relative weight is given by z„. A similar situation can be found 
in QED where the one-loop photon self energy is resummed by solving the Schwinger-Dyson equation. The infinite 
geometrical scries representation of the photon propagator can be interpreted as having infinite number of preons, 
corresponding to different number of particle-antiparticle pairs creation process. 



C. Preon condensate 



In order to see better the low energy preon dynamics let us consider a theory whose effective action, the generator 
functional of one-particle irreducible graphs of the field <j> is given by the space-time integral of the Lagrangian ([T])-®). 
The infinitesimal fluctuations of (j>(x) have relativistic dispersion relation with real mass m n , satisfying the equation 
An(m^) = 0. We restrict ourselves to the positive norm subspace where each preon has positive kinetic energy. 

Let us change slowly, adiabatically the effective potential V(<ft 4>) in such a manner that one or several m 2 approach 
zero and after that change sign. How does the vacuum change during this process? The particle number density, 
given by the propagator at zero separation in the leading order is the sum of the preon densities and the contribution 
of the light preons is obviously dominant when some of the preons become light. As soon as the preon mass squared 
becomes negative the vacuum energy is lowered by spontaneous preon-antipreon pairs creation until the repulsive 
interaction stabilizes the vacuum, a neutral coherent state. 

We look first for a static condensate and the possible time-dependent instabilities will be excluded by checking 
if small fluctuations have stable excitation spectrum next. The vacuum energy density is given by means of the 
expectation value of the energy-momentum tensor, (T 00 ) which is the sum of the tree-level and the loop contributions. 
The tree-level contribution, 

(T°°)tree = ^d Q <j> + d Q <jf J^-r - L + £ ho (11) 

contains a higher order derivative term, Eho- For static field configuration (T 00 ) tree — —L because space derivatives 
are left in this case only and they contribute to the potential energy. 

As long as preons with vanishing momentum build up the vacuum the condensate is space independent, homoge- 
neous. When preon states with nonvanishing momentum are macroscopically occupied then the interactions of preons 
in the vacuum can be taken into account on the tree level, by minimizing the volume integral of the static energy 
density. This minimization is highly nontrivial due to the competition of the space-derivative terms, (jy L(A)(j) and 
the local potential, Vfficj)). If L(— p 2 ) assumes its maximum at nonvanishing momentum then the configuration, 
minimizing the energy can be a space-dependent, periodic function. The simplest periodic vacuum configuration is a 
plane wave but oscillatory behavior can be observed in several directions, depending on the parameters of the theory 
according to the minimization of a 2-dimensional Euclidean field theory action, carried out in Ref. [43j which can be 
considered as the energy density of static configurations in 2+1 dimensional space-time. 



D. Goldstone preons 



It is natural to raise the question what happens when a continuous symmetry is broken spontaneously by the vacuum. 
Let us suppose for simplicity that the mass squared is non-negative for all preons and the vacuum is homogeneous 
and apply an infinitesimal symmetry transformation, <f> ~ ► 4> + <5</>j &4> — itTtfi, given in terms of a generator r* = r, 
which preserves the Lagrangian ([T]), 

V'(ft4>)ftr0 = 0. (12) 



6 



The derivative of this equation with respect to <p\ evaluated in the vacuum, 

O = M 2 t0 (13) 
is the usual Goldstone theorem for the mass spectrum. The dispersion relation of the Goldstone preons, 

det [L(p 2 )} = (14) 

M 2 =0 

the vanishing of the determinant within the null-space of M 2 contains a massless preon since L(0) =0. If p 2 = is 
a single pole within the null-space then we have a massless preon but there might be additional massive Goldstone 
preons, as well. If p 2 = is a higher order pole then the massless Goldstone preon is "confined" , it does not correspond 
to a particle-like asymptotic state and scattering experiments see no massless Goldstone mode. 

III. GAUGE MODEL 

The use of higher order covariant derivatives in gauge theories opens up new possibilities concerning space-time 
symmetries, namely the recovery of translation invariance of a certain family of inhomogeneous vacua and the spon- 
taneous breakdown of Lorentz symmetry [ll|. We shall discuss these issues in the framework of the model (TTJ by 
upgrading the global, semi-simple U(n) symmetry group of a local U(n) = SU(n) x U(l) symmetry. The Lagrangian 
of this Yang-Mills-Higgs model is L — L m + L g , where 

L g = -\t,F 2 -\ G 2 , 

L m = 4^L{-D 2 )^-V{^4>), (15) 

where the covariant derivative, — <9 M — igA^ — ieB^, is given by means of the gauge fields = A a ^r a and -B M , with 
[r a , r b ] = if abc r c and trr a T b = \8 ab . The field strength tensors are defined as F„ v = d^-dyA^-ig^, A v ] = F^r a 
and Gyu, = <9 M 5„ - d v B ll . 

Let us suppose that the vacuum of the model with global symmetry, g — e — is of the form (4>(x)) — uj(x)4>o, where 
uj(x) is a local U(n) transformation. The gauge transformation <p — > lo^4>, igA^ + ieB^ — > w(<9 M — igA^ — ieB^)uj^ 
brings this vacuum into the form (4>(x)) = 4>q, (igA^{x) + ieB^x)) — Lod^u:^ . A simple family of vacua which 
consists of a plane wave of an Abelian subgroup, ui{x) = e~ lkxT , it being a generator of U(n), yields {<t>{x)) = 4>q and 
(gA^(x) + eB^(x)) = fc^r. The dynamical mechanism of generating nonvanishing gauge field from a scalar condensate 
in the vacuum is provided by the higher order derivative terms in the kinetic energy of L m which contains a local 



potential for the gauge field [ll| with degenerate minima. The latter makes the interaction between gauge bosons 
attractive in the naive, "empty" gauge vacuum and the gauge boson density is stabilized by their repulsive contact 
interaction. 

Once a particle with nonvanishing spin forms a homogeneous condensate the invariance with respect to the Lorentz 
group is broken. This leads to the loss of symmetry under spatial rotations or Lorentz boosts when the condensate 
is space- or time-like, respectively. We explore below vacua with unbroken spatial rotational symmetry because this 
symmetry is expected to be present at any scale and space-like condensate violates time reversal invariance [lol ] . 

A. Mean-field vacuum 

The homogeneous, tree- level vacuum is based on the separation 4>{x) = (p + x(x), A^x) — A^ + a fJl (x) and 
B^(x) = B^ + bfj,(x) of the fields to a homogeneous condensate and quantum fluctuations. The basic building block 
of the kinetic energy of the scalar field, the argument of the function L can be written as the sum of the mean-field 
and fluctuation contributions, 



D 2 = □' - (aj, 6 M ) ( f g T) f a fe 7 ) + + bji (16) 



where □' = □ - {gA + eB) 2 - 2i{gA + eB)d, 



i a = -gr a {2id + i1) +2eB) -g 2 {T a ,A\, 

P = -e(2id + id +2eB + 2gA). (17) 
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The homogeneous, tree-level vacuum is found by minimizing the potential energy density 

U(A, B, = itr [{-ig[A^ A v ]f] - &L((gA + eB) 2 )4> + V($4>), (18) 

defined by omitting all space-time derivatives from the Lagrangian. We assume that the three-rotation remains a 
symmetry of the vacuum, hence the forms A^ = n^A, = n^B, with n 2 — 1 follow, reducing the potential energy 
density to 

U(A, B, = V{&$) - ftL((gA + eB) 2 )4>. (19) 

We shall consider the theory in static temporal gauge where the temporal components of the gauge field are 
static and behave as a three-dimensional adjoint scalar field with respect to the residual symmetry of this gauge, 
the space-dependent gauge transformations. Part of this symmetry is fixed by going into unitary gauge where these 
three-dimensional fields are diagonal, (gA + eB)ij = Sijkj. The notation fc' 1 = n^k will be used for the gauge field 
condensate, as well, k = gA + eB being an n x n matrix. The vacuum energy density can now be recast in the form 

n 

U(A,B,4>\4>)=V(ft4>)-Y / $L(k])^- (20) 

3=1 

The non-diagonal components of the gauge field are charged with respect to the remaining Abelian U(l) n gauge 
symmetry which is further restricted by taking the scalar condensate, <j>, real. 

We are interested in the case g A + eB therefore we need <j) 7^ 0. The gauge and the scalar field condensates may 
belong to different subspaces of the symmetry U(n). If m — Ka.ng(gA + eB) < n then we have the standard Higgs 
mechanism for n — m components of the Higgs field. We look for the interplay of the scalar and the gauge condensate 
therefore Rang(gA + eB) = n will be assumed in what follows. Minimization with respect to the gauge field gives 

L'(k 2 ) = (21) 
and L"(k 2 ) < as long as if>j ^ 0. The minimization in the scalar field yields the condition 

L(fc 2 ) = V'{ft4>). (22) 

Therefore, all kj components corresponding to <pj ^ must be at the maximum of L(p 2 ). The tree- level vacuum 
allows global gauge rotations as residual symmetries within the degenerate subspaces of kj . Different distribution of 
the magnitude of the scalar condensate among these degenerate subspaces yields inequivalent, classically degenerate 
vacua. 



B. Fluctuations 



The determination of the spectrum of infinitesimal fluctuations serves two goals, it establishes the stability of the 
static vacuum with respect time dependent perturbation and gives the particle content of the theory. 

It is advantageous to separate the real and imaginary part of the scalar field, x = Xl + *X2 and use the Fourier 
representation 



f(x) = / f(p)e-**, 



(23) 



with J — J <i 4 p/(27r) 4 for the fields. The calculation, sketched in Appendix \X\ gives the quadratic action 



g(2 



(Xi(-P), Xt{-P), aZ(-p), 6 M (-p)) 



(( L+ 



2V 



,,(^i®<t>i 0\ (\L d M hv \ [L d N"]+ 







r t [M^L d ] [-iM^La] 
5 i I {N»L d }+ -i[N»L d ]- 



Gg + G>, 



\ 


(xi{p)\ 




Xa(p) 






J \bu{p)J 



(24) 
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where G g and G m are given by Eqs. (|A4[) and (|A16|) . respectively, Ld = L((p + k) 2 ) — L(k 2 ) and 



.Y" = 



p 2 + 2pnk 
e(p^ + 2kn») 



(25) 



p 2 + 2pnk 

The use of real field components in space-time leads to the combinations f(p) + = \[f{p) + f(—p)] for real and 
f{p)~ = \\fip) — f(—p)] for imaginary matrix elements. 

We impose now static temporal gauge by writing = (a , na^ + ay), b^ — (6°, nfr^ + by), n = p/|p| in Fourier 
space. The quadratic action assumes the form 

/xi(p)\ 

S (2) = I (xi(-p),X2(-p),a a L (-p),b L (-p))K(p) g^j \+\J (a%(-p),b T (-p))G T (p) + (26) 



where 



tf(p) = 



V 



LfJ" 2y (, 
gr[r a JV£ 4i ] ± * fl r[7- tt JV£ e |] ± 



-ig[L d NT b 



61 ± 
± 



G 9 



-e[£ rf iV]+ 
-»e[iiJV]- 



(27) 



N = |p|/(p 2 + 2pnk) and denotes the contribution of the static modes, an expression similar to ((24)) except 
that the integration is over three momentum and the fields and the quadratic form is taken at vanishing frequency. 
The subscript ± of the complex matrix elements, acting on the non-Abelian gauge field is a reminder that the real 
and the imaginary components should be symmetrized or antisymmetrized, respectively. An important simplification 
comes from the fact that the higher order derivative terms in the kinetic energy, being constructed by means of the 
four- vectors p^ and do not act on the transverse gauge field components whose action comes entirely form L g , 



9" 



p^p" 



-5 ab p 2 



g 2 J deaf deb A c A e + 2igf abc A c n P 




(28) 



T J \ u ~P 

according to Eq. (|A4j) . 

The most important information in the quadratic action is the preon dispersion relation, given by the set of four 
vectors p^ = (w,p) where some matrix elements of the propagators K~ 1 (p) or G^, (jp) diverge. The inverse matrix 
element of a matrix is the ratio of a minor and the determinant. The former is a polynomial of finite order in p, hence 
the propagator can display singularity at the root of the fluctuation determinant, detK(p) or detGr(p)- 

The spectrum of the transverse gauge field components can easily be determined by noting that the action for these 
components comes from G g alone which is free of higher order derivatives. Therefore the gauge field condensate acts 
as a chemical potential and shifts the preon energies of the charged, non-diagonal gauge fields. It is advantageous to 
use the charge-diagonal basis, denned by [A, r a ] = q a r a because the quadratic form (|28|) is diagonal and the diagonal 
elements, (ui — gq a ) 2 — p 2 , give the massless spectrum, shifted by a chemical potential, oj = |p| + gq a ■ 

It is instructive to compare the spectrum with those, predicted with the conventional Higgs mechanism. The most 
obvious difference is that the transverse gauge field components remain massless, the term which is responsible for the 
gauge boson mass generation in the usual Higgs mechanism, (|A9j) . is absent from (|24]l due to Eq. ([2T]) . the possibility 
of accumulating gauge bosons in macroscopic quantity in the vacuum. Another remarkable feature of this generalized 
Higgs mechanism is that the preon spectrum breaks Lorentz invariance by a shift of the energy, a chemical potential 
only. 



C. Gap less scalar preon spectrum 



As the simplest use of the action (|26p we determine the gap in the preon spectrum. The gap is the energy at 
vanishing momentum hence we need the quadratic action for homogeneous, time dependent fields only. The static 
modes represent no dynamical degrees of freedom and can be neglected in (|26|) . leaving the quadratic action 

duo fL+ - 2V"6 1 (g 



2tt 

1 f dto 



-iL, 



iL d 



X2M 



+0 / ■^(fil(-u)M-<»)) 



LO 2 - 



g2 jdea jdeb J^c J^e 





2igf abc A c uj 



a b L (w) 



(29) 
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where the scalar and gauge fields do not mix. 

The action of the longitudinal gauge fields comes from G g and the remark, mentioned above about the transverse 
spectrum applies, giving the spectrum (uj — gq a ) 2 , displaying nonvanishing gap, ±gq a , for the non-diagonal, charged 
components. The scalar sector is gapless, as expected by Goldstone theorem, owing to the vanishing of the fluctuation 
determinant for u> = p = 0, 



det 



which is the product of the scalar preon masses. 



-2V>i<g)0i 








(30) 



IV. U(2) GAUGE THEORY 

It is instructive to consider the simplest nontrivial realization of the model discussed above, using the gauge group 
U(2), in more details and check the conditions of unitarity. We take the simplest choice of the scalar kinetic energy 
with a maximum, 

L(p 2 )=L - {p2 ~f f (31) 

K 

giving Ld(p) — — ((p + k) 2 — k 2 ) 2 / k 2 and kj — ±k. In case of identical signs, kj = (±k,±k) we have Abelian 
condensate, eB = ±«. When kj — (±k, =f«) then the gauge condensate is non- Abelian, gA a = 2/-c<5 a3 . There is no 
place for condensate containing both Abelian and non- Abelian components with such a simple kinetic energy and 
there are two, tree-level degenerate, inequivalent families of vacua. 



Let us assume that the gauge field condensate is non- Abelian, B^ 

'1 ' 



A. Non- Abelian condensate 

and 

/i n \ 

gA^ = 



-1 



K, 



(32) 



satisfying Eq. (12T1) . The scalar condensate (4>i, $2) — (sin 6*, cos 9)<fi, where real is chosen in such a manner that Eq. 
is satisfied. The vacua corresponding to different values of 9 are inequivalent and degenerate on the tree-level. 
The first term of the quadratic action, (|2T)|) . can be written as 



c(2) _ I f (h( _ ) of __^ ( K hhip) K hg {p)\ (h(p)\ 
~2j p [h{ Ph9{ P >> \K gh (p) K gg (p)) {g(p)J 

where h = (^x 1 , ^X 1 , ^X 2 , %X 2 ), 9 = (ah a 2 L , a\, b L ) and 



K, 



hg 



(2L+ - W"4> 2 sin 2 ( 


9 2iL- 




2V"4> 2 sin 26» 





\ 


-2iL- A 














-2V"(j> 2 sin 29 







- W"4> 2 cos 2 


9 -2iL~ d 


\ 





2iL- d 


2L 


\ ) 


/—gR~ cos 9 igR + 


cos 9 


-gR 


"sin^ -2eR~ 


' sin 0\ 




igR + cos 9 gR~ 


cos 9 


igR+ 


sin6» 2ieR + 


sin (9 




gR sin 9 — igR + sm 9 


-gR' 


'cos 9 2eRT 


cos 9 


- 


\igR + sm9 gR 


sin0 


—igR 


f cos 9 2ieR+ 


cos 9 ) 





K 



AK 2 +g 2 S + -Hun + ig 2 S~ cos29 



with 



AiujK — ig 2 S cos 29 


2egS+ sin 20 



R ± = 
S ± = 



'■+4k 2 + g 2 S+ 


-2iegS~sm29 



gh 



2eg5+sin26l 

2iegS~sm2e 

uj 2 +g 2 S+ -2egS+cos29 
-2egS+cos29 w 2 +4e 2 S*+ 



L d (p) 



± 



Ld(-P) 



(p 2 + 2lok) (p 2 - 2ujk) 
Ld(p) Ld(-P) 



± 



{p 2 + 2lok) 2 {p 2 - 2lok) 



(33) 



(34) 



(35) 
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The theory has unitary time evolution if all preon masses are real. The determinant of the quadratic form in Eq. 

dctK = -^cj 4 (w-2k) 2 (cj + 2k) 2 (w 2 -2cjk-p 2 ) 2 (cj 2 + 2wk-p 2 ) 2 

x [2K 2 4> 2 V"{2g 2 v 2 4> 2 - uj 4 + 2p 2 ui 2 - 4k 2 uj 2 - (p 2 ) 2 + 8e 2 p 2 <^ 2 cos 2 29) 



-ui 



4p 2 ^ 6 + 8k 2 lu 6 - 6(p 2 ) 2 ^ 4 - 16 K 2 p 2 ^ 4 - + 4(p 2 ) 3 ^ 2 + 8k 2 ( P 2 )W - (p 2 ) 4 ] (36) 



has complex frequency roots unfortunately. The source of trouble, the dependence on the scalar condensate, can be 
eliminated by imposing a further condition [Tl| . 

V" = 0, (37) 

when one finds 

detK = ^uj 4 {uj - 2k) 2 (uj + 2k) 2 (uj 2 - 2uik - p 2 ) 4 (w 2 + 2ujk - p 2 ) 4 , (38) 



and the preon spectrum is real and gapless, u) a . a > = an + & 'y 1 n 2 + p 2 , with a, a' = ±1. 

Let us now turn to the gauge boson propagation. The form (I28[) gives a massless transverse propagator with shifted 
energies for the charged gauge field components, 



G~^ v = ^ - ^ J I 1«-«-J'-p- ^ 1 . (39) 

written in the charge diagonal basis, a a = (a + , a~, a 3 ), with cP = (a 1 =p ia 2 )/V2 and <7 a = (2n/g, —2n/g, 0). To find 
the longitudinal propagator we first eliminate the scalar field from the quadratic action by means of its equation of 
motion. The elimination of x 1 and \ 2 generates the term 

5/ p (4(-,),«K-,))( ( "- 2 ' s)2 (lJ+ v)(:ig)' < 4o > 

with shifted energies for the charged fields. It also yields the contribution 

^[b L (-p)b L (p)+al(-p)al(p) . (41) 

f 

to the action. The quadratic action for a\ and fo^is diagonal when the field combinations 

g'b L + e'a 3 L 

ah 



ySg' 2 J- -' 2 ' 



'ff* + e' 
e'bL+g' 
y/g' 2 + ef 



* = ^'V^. (42) 



with g' = <? cos 2(9, e' = 2e are used, 



ui 2 [a L (-p)a L (p) + z L (~p)z L (p)}. (43) 



The theory, obtained by the fine tuning (|37[) is unitary in the positive norm, physical sector and has two charged 
scalar preons for each scalar field component, (j>j. At least one of these field components has gapless excitation 
spectrum according to Goldstone theorem. There are two, propagating massless transverse helicity preon modes for 
each gauge field a° and b^. The longitudinal gauge field components have a single, non-propagating and gapless preon 
mode as expected from the space-dependent residual symmetry of static temporal gauge. The energy spectrum of 
the charged components of the non-Abelian gauge field is shifted by the gauge field condensate as chemical potential. 
The absence of the parameter 9 of the vacuum in the eigenvalues of the propagator indicates that the vacua with 
different values of 9 remain degenerate on the one-loop level. 
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Abelian condensate 



Another family of vacua is where the Abelian gauge field develops vacuum expectation value, = and eB^ = n^K. 
The residual global U{2) can be used to rotate the scalar condensate into the form <j)j = (0,(f>), with real </>. The 
action. (1331) now contains the matrices: 



K, 



hg 



K aa = 

The fluctuation determinant reads 
16 



/ 2L+ 2iL~ \ 

-2iL- 2L+ 

2L+ - AV"{$) 2 2iL- 

' -2iL~ d 2L+J 

-gR- igR+ 

igR+ gR- 

gR- -2eR- 

-igR + 2ieR+ 



/uj 2 + g 2 S+ -ig 2 S~ 



^g 2 S- 





g 2 s + 
o 








-g 2 s^ 



-2egS+ u? + Ae 2 S 




(44) 



detK 



2ojk - 



2lok - \p\ 2 Y [2 K 2 (<P) 2 V"(2g 2 p 2 (0) 2 + 8e 2 p 2 (0j 2 - + 2p 2 u 2 - Ak 2 u 2 



-(p 2 ) 2 ) -uj» + 4p 2 u b + 8 K V - 6(p 2 ) 2 w 4 - 16k 2 p 2 uj 4 - 16 K 4 ^ 4 + 4(p 2 ) 3 W 2 + 8« 2 (p 2 ) 2 w 2 - (p 2 ) 4 ]. (45) 



It gives complex frequencies as in the case of non- Abelian condensate and we impose again the fine tuning condition 
(|57)) . yielding 



detK = -^(w 2 -2wk-p 2 ) 4 (cj 2 + 2o;k-p 2 ) 4 



(46) 



to establish unitary time evolution. The preon spectrum is the same as in the case of non- Abelian condensate, namely 
bj a (J i = <7K + (j 1 yj k 2 + p 2 . Any other L{p 2 ) = O (p 4 ) polynomial which is bounded from above and has maximum at 
k 2 reproduces the same fluctuation determinant up to a constant factor. 



The transverse gauge boson propagator, (|28|) is 



5 ab 
1 



The elimination of x 1 from the quadratic action by means of its equation of motion generates a term 



u 2 at(-p)a L (p). 



(47) 



(48) 



When \ 2 is eliminated then the contribution (|4Tj) is generated in the action. The full quadratic action for a\ and 
can be diagonalized by the help of the combinations (|42|) with 6 = and (|43|) is recovered. The spectrum is similar 
as in the vacuum with non- Abelian condensate except that the chemical potential which shifted there the energies is 
vanishing in this case. 



CONCLUSIONS 



Models with higher order derivatives are usually considered as effective theories, arising by the elimination of some 
heavy particles. A different motivation to consider such theories is followed in this work, namely the higher order 
derivative terms not only preserve renormalizability but act as a physical UV cutoff and opens the way to apply a 
wide family of new theories to realistic problems. 

The case of a spontaneously broken global gauge invariance is considered in this work where the higher order 
covariant derivatives in the scalar field kinetic energy generate a local potential for the gauge field which in turn 
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induces a nonvanishing expectation value of the gauge field. In order to preserve rotational invariance and unitarity 
only the temporal component of the gauge field can develop expectation value. 

In contrast to the usual Higgs mechanism the gauge field remains massless and the gauge field condensate acts as 
a chemical potential only. The mass term is removed by the stabilization of the vacuum with respect to the density 
of gauge bosons in the coherent state vacuum. In other words, the spontaneous breakdown of the symmetry with 
respect to Lorentz boosts consists of the identification of a preferred inertial system where preons have relativistic 
dispersion relation with a chemical potential and form a relativistic gas at finite density. 

The higher order derivative terms in the kinetic energy suppress the UV dynamics and generate several dispersion 
relations, preons, for a single field. Therefore these terms give a new class of unified models. For instance one may re- 
open the issue of massive gauge bosons without spontaneous symmetry breaking, Higgs mechanism and renormalizable 
models can be constructed where the price of massive vector bosons is the appearance of further preons, particle-like 
modes in the theory. 

Appendix A: Quadratic action 



The quadratic part of the action, defined by the Lagrangians (|15|) is calculated in this Appendix, following the 
scheme worked out for the Abelian case fill. 



1. L a 



Let us start with 



S a = -\ dxtr[(d^(A u + a v ) - d v (A^ + a M ) - ig[A„ + a„, A v + a u ]) 2 ] - - \ dx{d p b v - dyb^f ( A I ) 



which gives the quadratic part 

1 



Sf = I dx 



^tx[(d^a v - d v a^) ] + igtT[[a^, a v ]F^[A]] + 2igtv[{d^a v - d^a^A^, a„]\ 



+ytr[([4 1 ,o I/ ] + K,A v ]f] - l(d^K-d„b^ 2 
which can be written for the gauge field, given by the Fourier integral as 



4 2) = Jm-p)M-p))gT 



(A2) 



(A3) 



where 



GT = ^(pV - pV") + ( T a) (A4) 



and 



2 

Q{an)(bv) _ 9_f fcaf deb A c Ae(jl » n u + g ^ + l JL f^c A c [2{n . _ n „„ _ 

y 2 2 

2. Lm 

The scalar field kinetic energy contributes to the quadratic action by a number of terms which will be calculated 
by using the power series representation 

oo 

L{z) = J2 c n z n . (A6) 

The quadratic part in the gauge fields comes from 

^L(-D 2 )0 = Cnfihid - ieA - ieB) 2 ] n j>. (A7) 
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One kind of contribution is found by keeping the quadratic part of D 2 in the gauge fields, 

l E °» E • • • *H0) {ll f) (S) • • • (A8) 



which can be written as 



This term is vanishing due to Eq. (|2 1 1) . 

The other quadratic contribution of (| A7|) arises by keeping the linear part from two different covariant derivatives, 



71— 1 n 

+ b(~p)(3(~p)} -..(p + k) 2 --- [a a (p)a a (p) + b{p)p{p)] ■ ■ ■ k 2 4>, (A10) 



/E c «E E ftk 2 ---[a a (-p)a a (-p) 
J v n i=i e'=e+i 



where a a (p) = —g(r a p + {r a , k}) and j3 = —e(p + 2k), cf. Eq. (TTT1) . We rewrite this expression in a simpler manner 
by means of the functions X(p) = a a (p)a a (p) + b(p)(3(p) and fj m t = f{k 2 , (jp + k m ) 2 , k 2 ) as, 

E I jml4>\X- p j m X pm i(j>i (All) 



where 



/(a,6,c)=^c„^ £ J-W- 1 - 1 ^ ■ (A12) 

n 1=1 t'=i+l 



The sums can easily be calculated with the result 

b- 

(c-b)(c-a){b-a) 



L{c){b - a) + L{a){c - b) - L{b){c - a) 
f(a,b,c) = 7— -y- r (A13) 



for a ^ c which in the limit a, c — > k 2 becomes 



L(b)-L(k 2 ) 



by the help of Eq. <J2TJ). The O ((a + b) 2 ) contribution is therefore 



with 



(A15) 



3 M j$ (A16) 

and 

G%% {hv) = g 2 {r a p» + ^{r Q , fc}) jm (rV + rf {t\ k}) ml 
G%t V = 9e(T a p» + k}) jm (p v + 2kn v ) mt 

GfZ ] = eg(p^ + 2kn v )jm(r h p v + n v {r\ k}) ml 
GjL = e 2 (p^ +2kn fi ) jm (p u + 2kn u ) mi . (A17) 
The mixing between the scalar and the gauge field comes from 

L(-D 2 )4> = c ^~( d - W A - i9B) 2 ] n 4> (A18) 



whose linear part in the gauge field is 



^^(-□')-(aV + bp) •••(-□') 



n 1=1 

hence the O (xH a + &)) contribution of the scalar field kinetic term is 

n 



The summation over £ is trivial and we find 



CM) - E c »E a ^ 16 " 

n 1=1 

L{a) - L(b) 



yielding 



X l[L((p + k) 2 )Y p ~ Y p L(k 2 )}0 



with Ypjm = Xpj m /[(p + kj) 2 - k 2 n ] for (IA20|) . The contribution O ((a + 6)%), 

-^[K_ p L((p + fc) 2 )-i(/E 2 )y_ p ]x. 

can be found in a similar manner. 

Finally, the kinetic energy for the scalar field fluctuation assumes the form 

xH-p)L d (p) x (p) - \Wxi-p) + xH-pM^x(p) + x\p)^v"( 



with Ld(p) = L((p + k) 2 ) — L(k 2 ), after using Eq. 
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